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This paper makes a comparison between x-ray absorption (XAS) and resonant inelastic x-ray 
scattering (RIXS) in the rare earths. Atomic calculations are given for 2p — > 4/ and 2p — > 5d XAS. 
The latter calculation includes the contraction and expansion of the 5d orbitals resulting from the 
complete exchange interaction with the 4/ electrons. The radiative decay of the XAS final states 
is described for the situations where the core hole created in the absorption process is filled by a 
valence electron or by an electron from a shallower core level. RIXS spectra, 4/ n — + 3d4/ n+1 — » 4/ n , 
integrated over the outgoing photon energy (fluorescence yield) are compared with 3d — ► 4/ XAS. 
Sum rules related to XAS and RIXS and their applicability are discussed. 
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I. INTRODUCTION 

Experiments on rare earth systems have made an im- 
portant contribution to the development of magnetic x- 
ray dichroism and -resonant magnetic scattering. Both 
circularo and linearo magnetic x-ray dichroism have been 
observed for the first time in rare earths. The first satis- 
factory explanation of rcsaaant magnetic scattering was 
done for Holmium metal.EHa The electronic structure of 
rare earths is determined by the interaction between elec- 
trons in the localized 4/ orbitals and in the broad 5d 
band. Spectroscopies involving the 4/ shell can usually 
be successfully described by atomic multiplet theory. But 
whereas, e.g., magnetic x-ray dichroism at the M45 edges 
is rather well understood, the L23 spectra, dominated by 
the dipolar 2p — > 5d transitions, pose more problems. 
In the interpretation one has to take into account two 
effects. nn 

First, one observes pre-edge featureaim that are 
weak in the isotropic spectra but have a strong circu- 
lar dichroic signal. That these structures arise from 
quadrupolar transitions into the 4/ shell has been 
tablished by resonant elastic a^ir x-ray scattering): 
by the observation of a non-dinahir angular dependence 
of the circular-, dichroic XAS,Lr0 by resonant Raman 
spectroscopy^!!!] rand by partial deconvolution of the life 
time broadening .1111 

Second, the finite integrated intensity of the circular 
dichroism of the 2p — > 5d transitions results, not only 
from the polarization of the 5d electrons in the ground 
state, but also from a dependence of the 2p-5d radial 
matrix elements on the^direction of the 5c? moment rel- 
ative to that of the 4/.LE3 Therefore, band effects!] and 
the full (i/-Coulomb interactiorfjlij have to be included 
in the interpretation of the spectral line shape and the 
variations in the L23 circular dichroic branching ratios. 

A complication in the interpretation of RIXS is that 
the deexcitation cannot be simply decoupled from the 
absorption step as can be done when the excitation is 
far above threshold. Furthermore, the decay is different 
for transitions between two core levels or between the 
valence shell and a core level. 



For systems where a detailed description of the spectral 
line shape is complex, useful results can still be obtained 
from statistical methods. For XAS and x-ray scattering 
sum rules exist that relaia-the integrated intensities to 
ground state properties.Ejli3 These sum rules rely on a 
number of approximations. In general a constant radial 
matrix element is required. For sum rules that consider 
the spin-orbit manifolds separately one has to assume 
that the edges can be distinguished by the j-value of the 
core hole. For fluorescence yield no exact sum rules have 
been derived so far and the application of XAS sum rules 
requires numerical validationO In this paper we address 
the applicability of the sum rules. 

The paper is divided as follows. First, we show the 
separation of the geometric and electronic part for XAS 
and RIXS. Section III gives a derivation for the angular 
distributions for some common experimental situations. 
Section IV is devoted to the spectral functions and their 
sum rules. We describe the spectral line shape of the 
XAS spectra at the L23 edge. For the dipolar part we 
include the effects of contraction and expansion of the 
5d orbital by the df-Coulomb interaction. RIXS results 
are discussed for spectroscopies involving the 4/ shell. A 
comparison is made between the radiative decay in spec- 
troscopies where the intermediate state core hole is filled 
by a valence electron or by an electron from a shallower 
core level. We end with the conclusions. 



II. INTENSITIES 

The intershell transitions as a result of the absorp- 
tion or emission of a x-ray photon are described by 
Hint — ' A + A • p}. Expanding the vector poten- 



tial in plane waves gives 
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where is a normalization volume and ake annihilates 
a photon with momentum k and polarization vector e. 
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The absorption intensity is then given by Fermi's golden 
rule: 

I(kew) = \ \(™\H int \g)\ 2 6(Lu + E a - En)p u , (2) 
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with the density of oscillators given by p w — ^ . The 
different multipoles of the spectrum are obtained by ex- 
panding the plane wave in Bessel functions and spherical 
harmonics, i.e, 



e ikr = 5^[t]»Vt(fcr)k<*> .f (t) 



(3) 



with [a ... 6] = (2a + 1) ... (2b + 1). We use here the 
shorthands = k l C l (k) and for spherical tensors of 
rank one: k = k' 1 '. Note that k = 1. For fcr<d one 
has jt(fcr) = (fcr)*/[£]!!. With the use of the orthogonality 
relation one then obtains 



part, i.e. 



Iikeu) = ^-N u J2 B Q D tQ |<n|[e,k<*>]«. V*| 
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where cross terms between different tQ values have been 
omitted. The following factors have been defined: N u — 



n%r^> "Q ~ (f) 2 gives the relative tran- 



sition probability, where for the dipolar and quadrupo- 
lar contributions one has B\JB\ = -|(-§) 2 ; and D t q = 
[tQ]k t /{[t}\\btQ{k)B Q ). The different multipole spectra 
are given by 



I 4 «» = -£ J-pQ*(gnng) 
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(4) with £ n = uj + E„~E n + i^ where T is the lifetime broad- 



where the tensor couplings are defined in the Appendix. 
It is convenient to define the operators 



bt Q (k) 



{ [p,rW]« + (-l)«[rWp]« }, 



Q\-(t) „iQ 



(5) 
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where the factors btQ,(k) will be chosen in such a way 
that nicely defined operators are obtained, see Table Q. 
For tQ = 01 (electric dipole) and boi(k) — \/3 we have 

(n|V 01 |.g> = -^-(n|p|g) = ±{ n \[g-,r]\g) = (n\r\g). 
V6imui nui Zra 

By using the definition for the outer product and bn (k) = 

k, we find that V 11 is equal to ^ifp,^ 1 = f ^L, 

with a the fine-structure constant and ao the Bohr ra- 
dius. V 11 forms together with the gsS term the mag- 
netic dipole operator. Magnetic dipole transitions are 
about (a/2) 2 , i.e., five orders of magnitude smaller than 
electric dipole transitions of the same wavelength. Fur- 
thermore, we have for the electric quadrupole operator, 

V 12 =r (2) j that &i2 ( fc ) = 

The multipole expansion enables us to separate the 
absorption intensity into a geometric and an electronic 

TABLE I. Relevant constants for electric dipole (iQ=01), 
magnetic dipole (11), and electric quadrupole (12) transitions. 
The factors btQ(k) appear in the definition of the operators 
y t( 2 £,2^ gj ve relative probabilities of dipolar and 

quadrupolar transitions. The factors D t Q are used in the 
definition of the angular dependence T t( ^ z . 



Q V tc 



b tQ (k) 



D 



tQ 



el. dip. 
magn. dip. 
el. quad. 



a a j 



:',() 



3 

3^2 

Mir -sV2 



3 
1 

3 

1 lk\2 



cning of the XAS final states; the different combinations 
of the matrix elements are defined as 



V^(abcd) = [<a|(Vt)"2|6),(c|V*Q|d>]W , 



(8) 



where the normalization constants hq z are defined in the 
Appendix. The angular distribution is given by 



T^ z (ke) = [z]D 2 tQl 



The RIXS intensity is proportional to 



(9) 
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In the remainder of the paper we consider the different 
multipole transitions separately; to reduce the number 
of indices we remove the tQ of the transitions. In a way 
similar to that for XAS we find for the RIXS cross section 
with Q-polar excitation followed by a Q'-polar deexcita- 
tion 

/(kew.kVu/) = N u i Bqi N u Bq 



T zz ' r (ke, kV) ■ l zz ' r (uj, uj'). (11) 



The electronic part has now become a tensor product of 
absorption and emission 
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i**W) = mJEtp: 
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where £ 



E n 



E 



j with 7 the lifetime 
broadening of the x-ray inelastic scattering final states. 
The angular dependence is given by 



T zzr (ke,k'e') = [T 2 '(k'e'),T 2 (ke)] r ? 



(13) 
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III. ANGULAR DEPENDENCE 

The tensors for the angular dependence are chosen in 
such a way that T°=l for tQ=01, 11,12. For electric 
dipole transitions T z is k independent, and the C = 
component is, with respect to the Z-axis of our system, 

T z (e) =3[z](-l) 1+ *n Qz n llz U llz (e*,e,Z), (14) 

where the bipolar spherical harmonicsE3 are given by 
U xyz (a 7 b,c) = (-)y +z n-^ z [a^b<-y^ z -c^\ The bipolar 
spherical harmonics with xyz relevant for dipolar transi- 
tions are given in Table ||; expressions for higher values of 
xyz are given by Thole and Van der Laan.EU For general 
multipole transitions it is more convenient to recouple 
the angular dependence 

'lxl 

T z (ke) = (-l) 1 + t + z D*Y\xyz]{ t y t 



xy 



x [[e%ef,[k« k«]Tn Q . 



(15) 



It is straightforward to show that for left and right cir- 
cularly polarized light one has 

1/1",* , II i / i\i+mr,» _ -ia=\ J a xRn x C (k) TO=0 

I I m=l 

with a x = 1, |i, —\ for x = 0, 1, 2. After that one finds 
with the use of the theorem to couple two spherical har- 
monics of the same vector (see Appendix) that the angu- 
lar dependence is a function of C z (k)Jl 

i(T*(ke + ) + (-ir+™T*(ke_))=( TzC2(k) m= ° (16) 
I m=l 



where the coefficients 



%=(- 



I x I 

u,\ t y t } n Qz n llx n tty n xyz 
Q z Q 



are given in Table [II 



IV. SPECTRA AND SUM RULES 

A. XAS into the 4/ shell 

Here we consider two absorption edges often studied: 
the L23-edges, corresponding to transitions from the 2p 
orbital into the valence shell and the M45-edges that are 

TABLE II. Bipolar spherical harmonic U xyz (a, b, c) rele- 
vant for the angular dependence of dipolar transitions. 



dominated by the dipolar 3d — > 4/ transitions. The 
electric dipolar and quadrupolar transition operators are 
given by the Wigner-Eckart theorem: 



7 Aer 



\v q \ g )=Pci^2(-y- x ( \ Q q ( ; ; ) <ni<l«- 



where l\ creates an electron in shell I with orbital com- 

Acr 

ponent A and spin a. The reduced matrix element is 
given by 

PcliQ) = {-l) l ni Q c[lc}i J drR nil (r)rQR ncC {r). (17) 

For rare earths the crystal fields on the 4/ electrons 
are often negligibly small and one can assume spheri- 
cal symmetry (SO3). For a magnetic system the sym- 
metry is lowered to SO2 and the J-values branch into 
Mj = — J, . . . , J. Here we take the magnetic axis along 
the Z-axis and the ground state to be Mj — J. To ob- 
tain a non-zero intensity the total transition operator has 
to be totally symmetric. This means we have to consider 
1$ (uj) for XAS and lf r (u, w') for RLXS. I§ is a combina- 
tion of the matrix elements I q (gnng) = (p\VJ \n) (n\V a \ g). 
The combinations for Q = 1,2 are given in Table IV . 
For dipolar transitions we have the well-known spectra: 
isotropic (I®), circular dichroic (/q), and linear dichroic 

(/ ° 2) - nn 

Figures |l|-|j give the XAS spectra for transitions into 

the 4/-shell at the L2.3 and edges. Calculations, 

were done in the atomic limit using Cowan's programs.Ej 
The Hamiltonian includes the Coulomb interactions in 
the 4/ shell and those between the 4/ shell and the core 
hole and the spin-orbit coupling. Parameters were ob- 
tained in the Hartree-Fock limit and the values for the 
Coulomb interaction were scaled down to 80 % to ac- 
count for screening effects. The zero of the energy scale 
corresponds to the energy of the lowest eigenstate in a 
spin-orbit manifold. Figure [l] gives the isotropic 2p — * 4/ 
XAS spectrum. The used Lorentzian of 2 eV is smaller 
than the expected 2p lifetime broadening. pThis was done 
since recent experiments by Loeffen et alii demonstrate 
the possibility of partial deconvolution of the lifetime 
broadening. By deconvoluting high-quality XAS data 
by a Lorentzian with a width of 3 eV, they obtained 
a clear separation of dipolar and quadrupolar features. 
Also the spectral line shape of the pre-edge structures was 
in agreement with multiplet calculations. This method 
seems to be well suited to 

TABLE III. Coeficients T z defined in Eqn.@ for Q= 1,2 
and z — 0, . . . , 2Q. 
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FIG. 1. Isotropic quadrupolar 2p — > 4/ XAS spectra as a 
function of excitation energy uj for trivalent rare-earth ions. 
The left and right panel show the L3 and L2 edge, respectively. 
Spectra for the L2 edge have been scaled by a factor two. 

study the relative contributions of quadrupolar and dipo- 
lar transitions. The spectrum obtained after taking the 
difference between left and right circularly polarized light 

TABLE IV. The spectra I§ expressed in I, for Q=l,2. 



Ig=h+I +I-1 

Q=l (dipolar) Jq 1 = h - J_i 

I% = h-2I + I-i 

Io=l2+h + Io+I-i + I-2 

Io = I2 + — h^- 1 ~ 

Q—2 (quadrupolar) Iq — I 2 — \h — Io — §J-i + I-2 
1$ =h- 2/i + 2/-1 - J_ 2 
Io = h ~ 4/i + 6/ - 47_i + I-2 
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FIG. 2. Quadrupolar 2p — > 4/ XAS spectra as a function 
of excitation energy u> for trivalent rare-earth ions. The spec- 
trum for left minus right circularly polarized light is a angular 
dependent combination of I 1 (solid line) and J 3 (dotted line), 
see text. The left and right panel show the L3 and L2 edge, 
respectively. Spectra for the L2 edge have been scaled by a 
factor two. Furthermore, some of the spectra are scaled rela- 
tive to the spectra of Fig. |l[ The scaling also applies to the 
right panel. 

is a combination of I§ with odd z. Using Eqn. ( |l6| ) we 
obtaintl 

\{I(ke+uj) - /(ke-w)} = ^N^BH-C^lHu) 

+C 3 (k)/ 3 H}, (18) 

where C§(k) are Legendre polynomials of order z, i.e., 
P z (cos6?) with 8 the angle between the direction of the 
light and the magnetic axis. The spectra Iq(u>) and Iq (uj) 
are given in Fig. p|. 
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FIG. 3. Isotropic dipolar XAS spectra (4/™ -> 3d4/ n+ \ 
dotted line) and fluorescence spectra for isotropic incoming 
and outgoing light (4/ n 3d4/ n+1 -» 4/™, solid line) as a 
function of excitation energy w for trivalent rare-earth ions. 
XAS and fluorescence have been normalized to the intensity 
integrated over both edges; the scaling between spectra for 
different ions is arbitrary. The left and right panel show the 
M5 and M4, edge, respectively. 

Figures ^ and ^ show, respectively, the isotropic and 
circular dichroic M^g XAS spectra, |X=0.8 eV). These 
spectra have been published earlierc3~Ej and are included 
for comparison with the fluorescence yield spectra. 

The behaviour of the intensities can be -tedated to 
ground state expectation values by sum rulesJia£3c3 For 
I = c + Q one can reduce the expressions tcOeS 

F {j) = (<L>F {u>) = ^X2 M y(j)N xyz (< yz ) (19) 

3 L J X y 

with M v (j)=c + l,c, c,-c for jy=j + 0, j~0, j~l 
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FIG. 4. Same as Fig. g but now for left minus right 
circularly polarized incoming light. Some of the spectra are 
scaled relative to the spectra of Fig. |§|. The scaling applies 
also to the right panel. 

(j* = I ± s) and N xyz =l,-fa,^ for xyz = zOz; z - 
1,X,jsiz+1, 1, z. The coupled tensor operators are defined 
asOEl 

_ ! \l-x'+s-<r ( 1 x l\ ( s y s\ 

W C - 2^ y > I -A' £ A ) \ -a' V a ) 

x (-) 2 - C ( ^ \ I ) h'A* «R- V \ (20) 

with s=i. The operators are spin independent and de- 
pendent for y = and 1, respectively. These hole cou- 
pled tensor operators are related to the electron oper- 
ators defined by Carra et alxQ via O xyz =2[l]8 x ,o<>y.o — 
{-l) z ri x r sy WQ yz with r/ x = 2 *(2i-xy. ■ Tbe advantage of 
using normalized operators over Judd's operatorsE3 is 
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FIG. 5. The upper panel shows the ground state expec- 
tation values of w z0z for trivalent rare-earth ions for z= 1 
(solid line), 2 (long-dashed), 3 (dashed), and 4 (dotted). 
The lower panel shows the ground state expectation value 



of = ft, 



+ ij±^io*+ 1 ' 1 '* (solid line). The dashed 

and long-dashed lines show rih[Io(j + ) — Sl £~Io(j ")] for the 
L23 and M45 edges, respectively. This quantity is equal to O z 
according to the sum rules, see text. The thin and thick lines 
correspond to z=0 and 1, respectively. 

that the square roots are removed from the expressions. 
They are normalized in such a way that the expectation 
values are unity for a ground state with one hole and 
Mj = J = I + i (i.e., for rare earths the ground state 
of Yb 3+ , f 13 ( 2 Fr;Mj = |)). Some typical examples are 
the number of holes: rih 
magnetic moments: L z 

tively, and the spin-orbit coupling: L S = —Iswq W . Note 
that the relation for O xyz does not hold for the magnetic 
dipole operator T z — ^rp^o 11 due to its definition in 
real instead of angular momentum space. 

By summing Eqn. (19) over both spin-orbit split edges, 
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swq 11 , respec- 



Io(3 + ) + Io(D 



(21) 



we obtain the sum rule that relates the totaLintegrated 
intensity to the spin-independent operators.li3 (I indi- 
cates that the spectra are normalized to the isotropic 
spectrum). For z — 1 we have the "L 2 " sum rule. In 
general, it relates the integrated intensity of I§(u>) to 
spin-independent operators that give the 2 z -polar mo- 
ment in the electron distribution. For z = it simply 
says that the integrated intensity over both edges is pro- 
portional to the number of holes. For higher values of 
z the ground state expectation values of w z0z show an 
oscillatory behaviour along the rare earth series, see Fig. 
H This behaviour is alreadv-well described by assuming 
a Hund's rule ground state.ea 



A weighted substraction of both edges can be related 
to the spin-dependent operators, 
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For z = this expression forms the basis of the theory of 
branching ratios It relates the integrated intensities of 
the spin-orbit split edges of the isotropic spectrum to the 
ground state expectation value of the spin-orbit operator 
(note thajUto- 1 ' 1 ' = 0). For z = 1 it gives the "S z -T z " 
sum rule.E£l 

The assumption in deriving this equation is that a cer- 
tain spin-orbit split edge can be described by the core 
hole having the corresponding j-value. This approxima- 
tion is valid for late rare earths (and also late transition 
metals). For early rare earths deviations are found for the 
spin-dependent sum rule, see Fig. |5|. Although for these 
systems a smaller spin-orbit coupling is found, this is not 
the dominant effect that causes the stronger mixing of 
the two edges. More important are the LSJ values that 
are reached by the absorption. For late rare earths the 
ground state has maximum LSJ values. As a result of 
the dipole selection rules (AS = 0, AL = 0, ±1) also high 
LSJ values for XAS final states are found. These states 
are predominantly found in the j + -edge. A qualitative 
argument for this goes as follows. Within LS'-coupling 
(i.e., no spin-orbit coupling) the cl n+1 states with high 
LS values have relatively low energies (where c indicates 
a core hole). After switching on the spin-orbit coupling 
the LS states with low energy predominantly go into the 
j + -edge and those with high energy into the j _ -edge. 
Hence, for the states with high LSJ values there is only 
small mixing between the two edges. The most extreme 
examples are the states with maximum LSJ values that 
only occur in the j + -edge. The limited presence of high 
LSJ character in the j~ edge directly explains why the 
the L2 and M4 edge have little intensity for late rare 
earths. 

The situation is different for early rare earths. First, 
the total J of the ground state is given by \L — S\. 
Second, for early rare earths the maximum spin of the 
XAS final states is .higher than that of the ground state 
(5™ ax = S™ 3 * + 1)B Since in XAS AS = these states 
have little weight (as a result of the mixing by the spin- 
orbit coupling they obtain a finite intensity). For early 
rare earths excitations therefore occur at intermediate 
LSJ values. This can also be seen in the M45 XAS spec- 
tra where for early rare earths significant intensity only 
occurs at 5-10 eV above the absorption edge (which con- 
sists of the "dipole-forbidden" maximum spin states) , see 
Fig. |^. The stronger mixing between the edges is directly 
apparent from strong intensity at both edges. This mix- 
ing has as a result that for early rare earths deviations 
for the spin-dependent sum rule are found, although for 
z = 1 the trends are well predicted. 
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B. 2p ~ » 5d XAS and elastic scattering 

We now discuss the behaviour of XAS and resonant 
elastic scattering at the 2p — > 5d absorption edge. For 
clastic scattering, it is more convenient to couple the ma- 
trix elements of the excitation to those of the deexcita- 
tion instead of coupling the excitation and deexcitation 
to their complex conjugates, as was done for inelastic 
scattering, see Eqn. (Jll|). One then obtains for dipolar 
transitions 



Y,^T*(e,e')-r(gnng)f 



f 



where the angular dependence T z (e, e') is given by Eqn. 
(fll|) under the replacement of e* by e'*. For XAS, the 
z = 1 spectrum is selected by taking the difference be- 
tween left and right circularly polarized light. In elastic 
scattering one adopts a a — > tt scattering geometry, where 
a and n denote linear polarization perpendicular to and 
in the scattering plane, respectively. From the expres- 
sions for U llz (e a , e„, Z), see Table Q, one sees that the 
z = 1 contribution is obtained if Z, i.e., the magnetic 
axis, is in the scattering plane. Thus the a — > n spec- 
trum can be approximately considered as the square of 
the circular dichroic XAS spectrum. 

The difficulty in the interpretation of the 2p — > 5d cir- 
cular dichroism is that the radial matrix elements are 
strongly affected by the presence of 4/ electrons. Since 
the radial extent of the 5d orbital depends on the di- 
rection of the 5d moment with respect to that of the 
4/ (and therefore on energy) deviations from the XAS 
sum rules should be found. This was indeed observed 
experimentally where the sign of the signal could not be 
reconciled^ with the expected direction of the magnetic 
moment .ElrE3 The change in radial matrix clement was 
first considered by band structure calculations.!] How- 
ever, in LDA the circular dichroic branching ratio is 
always —1, since only the bd-Af spin interactions are 
considered. Recently, a model has been developed that 
includes the complete d/j-Coulomb interaction, i.e., also 
the orbital contributions]!.] The assumption was that the 
dichroic signal was mainly determined by the change in 
the radial matrix elements and that the 5c? polarization 
in the ground has a relatively small effect, i.e., one can 
assume a bd° ground state. For constant radial matrix 
elements, the integrated intensity of the circular dichro- 
ism is then zero.t3'c3 However, this does not mean a zero 
signal. As a result of the interaction with the 4/ elec- 
trons the circular dichroic spectrum contains equal pos- 
itive and negative parts. Note, that this implies a finite 
elastic scattering amplitude even in the absence of a bd 
polarization in the ground state. Due to the large 2p 
lifetime broadening the spectrum for g-polarized light is 
mainly determined by its first moment, (g\VJ HdfV q \g) . 

Within a framework with fixed atomic orbitals one can 
create a bd orbital with a different radial extent by mixing 
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FIG. 6. Dipolar circular dichroic 2p — * 5d XAS spectra as 
a function of excitation energy uj for trivalent rare-earth ions. 
The left and right panel show the L3 and L2 edge, respectively. 
All spectra have the same scale. 

it with other n<i orbitals, i.e., |5rf) = a\5d)+^2 , g a n \nd). 
Using first order perturbation theory one findsEJ that the 

zeroth moment I§ of the absorption spectrum into 5<i is 

{1)7 

proportional to the first moment Iq of the absorption 
spectrum into 5d: 



n^5 



n^5 



Y,L kx (df)M v (j)N xl 



(23) 



where A n d is the energy difference between the bd and a 
nd orbital, and 



7 



L kx (df) = (-) 



If] 



dkf 



c n-dx 



f f k 
d d x 



(24) 



The F k direct Coulomb terms have been omitted since 
they do not contribute to the circular dichroic first mo- 
ment. Note that the ground-state expectation values in 
Eqn. ( p3| ) are those of the 4/ and not the 5d shell. The re- 
sult looks very similar to the XAS sum rules of Eqn. ( 19 ) . 



In particular, for the k=l terms one has L^ x {df) = 1/35 
and, except for a scaling, the same coefficients as for XAS 
are found. 

Figure || shows numerical calculations of the 2p — > 
5d circular dichroism at the L23 edges using final 
state configuration interaction between 2p4f n+1 5d 1 and 
2pAf n+1 nd 1 . To reduce the size of the calculations only 
the 6d orbital is included. To account for the other nd 
orbitals the Coulomb interactions between the 5d and 
the 6c? shell are multiplied by four. The energy difference 
between the 5d and 6d orbitals has been taken 15 eV; 
only the 5d region is shown in the figure. The broaden- 
ing is 6 eV which accounts for band effects and lifetime 
broadening. As a result of computer limitations we were 
unable to obtain results for Gd 3+ . One would expect the 
spectrum to be given by a single negative peak at the L3 
edge and a positive peak with equal intensity at the L2 
edge. 

Despite the crudeness of the model one finds-qualita- 
tive agreement between theory and experiment £2rE3 For 
late rare earths the signal at the L3 edge is larger than 
that at the L2 edge, whereas for early rare earths the 
L 2 is the larger. The sign of the Iq quadrupolar contri- 
bution to the absorption spectrum is opposite to that of 
the dipolar, except for the L 2 edge for 4/ 9 to 4/ 13 where 
the same sign is found. This difference is an effect of the 
G 3 and G 5 terms in Eqn (p3[). Also the relative signs of 
the edges are well predicted. For early rare earths both 
edges have the same sign, whereas for late rare earths an 
opposite sign is found. 

One also observes that, for early rare earths and for 
the L2 edges in late rare earths, the spectrum is more 
derivative-like than for the L3 edge in late . rare earths. 
This effect is also seen experimentally.EJO We do not 
have a satisfactory qualitative explanation for this effect. 
We also note that for the L 2 edge the signal occurs at 
lower energies with respect to the lowest final state. This 
effect, combined with the larger lifetime broadening of 
the J 2 -edge, might explain why the quadrupole contri- 
butions are not always easily observed in the L2 edge. 

For finite bd polarizations one can approximate the to- 
tal intensity by I§ + Iq, where I§ is given by Eqn. (19). 



Note that the I§ and I§ are given by 5d and 4/ ground- 



state expectation values, respectively. Experiment shows 
that Iq dominates. The trends in \Iq | are similar to those 
of I Jo I, but their signs are usually opposite.El A finite 5d 
polarization therefore decreases the integrated intensity 
(in principle it could even change sign)^ but has rela- 
tively little effect on the branching ratio. 



C. Resonant inelastic x-ray scattering 

For x-ray inelastic scattering sum rules one has to 
integrate over the incoming and outgoing photon ener- 
gies. Let us first consider the integration over the trans- 
ferred energy, Aw = w — u', which leads to a term de- 
scribing the decay of an intermediate state, (n'n) — 

J2f Ifi ( n> f f n )- Note the presence of cross terms as a 
result of interference between intermediate states lead- 
ing to the same final state. For the decay one has to 
distinguish two situations. 

First, the core hole is filled by an electron from a shal- 
lower core level, leading to a Raman process described 
by l n — > cl n+1 — > c'l n+1 . If both the intermediate and 
final states are split by the spin-orbit coupling the spec- 
trum has four clearly separated manifolds. For a given 
manifold we findcS 



I* c ;(nn';jj')=Bf f (c,Q,c')J2H 



j—m' [ J % J 

—m! Cj m 



x(n'\c jm ,c] m \n)nj, , (25) 



with the coefficient given by 



B(Ac,Q,c>) = (-y + 3' +c+c 'P? c ,[jf 



J 3 J' Q' 



, Q Q j' i n ., n -i 

' 3 3 z' 1 ]Z Q' 2 



(26) 



The decay of an intermediate state (n = n') is therefore 
determined by an expectation value of the polarization 
of the core hole for that state. For, e.g., z' = 1 this 
expectation value is (j z )- 

Let us considep-a specific example. A typical resonant 
Raman processcHl3 is given by /=4/, c=2p, and c'=3d. 
For simplicity we take the incoming light isotropic. If the 
system is magnetic the absorption process still creates 
a polarized core hole as a result of the polarization of 
the valence shell. If the polarization is not detected the 
outgoing light has an isotropic and linear dichroic part, 
see Eqn. (|l6|), 

k'e+w') + 7(w, k'e-w')} 
= 2 -^N„,N u BlBl{ir(u,u') + icg&ifVX)}- 

Figure gives the two different spectra for Ho 3+ at three 
different incoming photon energies. If we now take, for 
a certain absorption energy, the difference of two spectra 
at different detection angles with respect to the mag- 
netic axis, we obtain a signal that is proportional to 
Iq 22 (uj, uj'). In that case the integration over to' of a spec- 
trum at a given absorption energy is finite at the L3-edge, 
but zero at the L 2 -edge, see Fig. [?|. This is a direct re- 
sult of the fact that the j ' = | level has no quadrupolar 
moment (a finite intensity can still occur as a result of 
the presence 
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FIG. 7. RIXS spectra (4/ n -> 2p4/ n+1 -> 3d4/ n+1 ) as a 
function of transferred energy Acj = u — u>'. The zero of 
the energy scale corresponds to the energy of the lowest final 
state of a certain edge. Spectra are given at three excitation 
energies, from top to bottom: at the Lz absorption edge, 2.5 
eV above the Lz absorption edge, at the L2 absorption edge, 
see Fig. |l[ Only the M5 edge is shown. The incoming light 
is isotropic; the outgoing light is isotropic (Jq (w, lu'), solid 
line) and linear dichroic (Iq 22 (u>, a/), dashed line). The dotted 
line gives the intensity of Iq 22 (us,u') integrated along the Auj 
axis. 

of j = I character in the Li edge due to the mixing by 
the core- valence Coulomb interactions). 

For a measurement with isotropic outgoing light the 
expectation value is proportional to #n,n'> i-e., the decay 
is constant and all interference terms cancel. WfecH-,a 
summation over both j'-edges is made one obtainal3'E2l 

/jj'=°(nn') = 5>o°K'; jj') = ?ff5 n<n , . (27) 

Therefore, for isotropic outgoing light an integration over 
the transferred energy u> — uj' leads to a spectrum as a 
function of excitation energy uj which is proportional to 

xas.be! 

The second case is where the deexcitation involves the 
valence shell, i.e., l n — > cl n+1 — > l n . Figures || and |] 
show the spectra integrated over the energy of the outgo- 
ing photons (fluorescence yield) for isotropic and circular 
dichroic incoming light, respectively; the outgoing light 
has been taken isotropic. We observe significant differ- 
ences with the XAS spectra. Therefore, even for isotropic 
outgoing light, the decay cannot be decoupled from the 
excitation. In contrast to Eqn. (|25|), where the operators 
working on the closed shell could be removed, one obtains 
here that the decay is proportional to an intermediate- 
state expectation value of a tworparticle operator. For 
isotropic outgoing light one hasia 
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x c 1 {cj',l\')c 1 {cj,l\)(n'\l x<J cy a d{, (7 ,cl (<T \n)}, 

where c 1 (07, IX) are the Slater-Condon parametersS The 
matrix elements are, except for an offset and the radial in- 
tegrals, equivalent to those of the G x cl Coulomb exchange. 
Note that the a G\" term also assumes negative values 
and that the lowest possible value for I§ =0 is zero. The 
G x cl Coulomb term is also for a large part responsible for 
the position of the eigenstates in a spin-orbit manifold. 
This directly explains the trend in the fluorescence spec- 
tra, that the states closer to the absorption threshold in 
a certain edge have in general a smaller decay compared 
to those at the high-energy side of an edge. 

Let us now consider the sum rules for the intensities 
integrated along the transferred and excitation energy. 
For deep-lying core levels with a large lifetime broaden- 
ing, such as tha,2» shell, one can use the fast-collision 
approximationOila This implies replacing the interme- 
diate state energy denominator by £ n = uj + E g — E n + 
i^, where E n is an average energy. For the situation 
where, the. emission involves two core levels one can then 
deriveM 

lf r {j,3')=Jdwj^'lf r {^') 

= Bjj, (c, Q,c')J2 C'7 rzz ' (c, Q, I) K yr > (28) 

xy 

where the coefficients B and C are given in Ref. ^8|. For 
isotropic outgoing light we have Cj VZz0 — ^M y (j)N xyz , 
which means that the same coefficients are found as for 
XAS. In this limit the fast-collision approximation is not 
necessary since all interference effects cancel. 

When the deexcitation involves the valence states the 
situation is more complex. Sum rules would involve two- 
particle valence shell expectation values which are dif- 
ficult to evaluate. Recently, an investigation has been 
made into the applicability of XAS sum rules for spectra 
obtained with fluorescence yieldJij Although fluorescence 
yield is in principle not equal to XAS the conditions for 
integrated intensities are less stringent. Here one does 
not require that every final state has a constant decay 
but that the total decay of the excited intermediate states 
does not have a strong polarization dependence. 

The state after excitation with g-polarized light can 
be written as \v q ) = J2 n (<?) I n ) ! this state can be 
normalized to unity by using the coefficients a n (q) — 
a n (q)/ yj (v q \v q ). The integrated intensity of the fluores- 
cence spectrum excited with q-polarizcd light can then 
be written asEj 

^ fluor = (v q \v q )(V r ) q = lf AS (V r ) q , (29) 

i.e, the integrated intensity is given by the XAS intensity 
multiplied by a term that describes the radiative decay. 
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If the polarization of the outgoing light is not measured 
this latter term is given by 

z' even nn' I — ^-p — — 1 -|- 1 

This clearly shows that the total decay for g-polarized 
light {V T ) q is a weighted average of the decays of the 
intermediate states I§ (nn 1 ). Proportionality of the inte- 
grated intensity of the circular dichroic fluorescence spec- 
tra, i.e., (£ ff qI^ or )/(E q ^ u ° r ), with (w 101 ) = (L z )/l is 
obtained if (V r ) q is not strongly dependent on the polar- 
ization. Numerical evaluation of (V r ) q shows that this 
situation is found for early rare earths, but that there as, 
a strong polarization dependence for late rare earths S3 
This can be understood as follows. As was shown above, 
the decay of the intermediate states is proportional to a 
"Gj ; "-like term. Maximum variations in decay can be ex- 
pected for "pure" LS-like intermediate states. The effect 
of mixing by the spin-orbit coupling is small for the high 
LSJ states that are reached in absorption in late rare 
earths, leading to a strong polarization dependence. The 
stronger mixing of the intermediate LSJ-states that are 
reached in early rare earths decreases variations in (V T ) q . 
It is remarkable that the mechanism that causes devia- 
tions for the "S z " XAS sum rule improves the agreement 
between L z values obtained by fluorescence yield and the 
ground-state expectation values. 

It must be noted that a more efficient way to re- 
move variations in (V r ) q is to create a \v q ) which has 
not a strong LiS-like character. This is found for many 
transition-metal systems where a strong crystal field 
quenches the orbital moment. 

V. CONCLUSION 

In conclusion, a comparison has been made between 
XAS and resonant inelastic x-ray scattering. Several as- 
pects deserve further experimental investigation. The 
isotropic L23 spectra should be reexamined since, as 
was shown by Loeffen et aZ.,EJ more detailed information 
can be obtained by partial deconvolution of the lifetime 
broadening. This should enable a direct comparison of 
the relative size of the dipolar and quadrupolar contri- 
butions. Also a careful determination of the intensities 
of the two spin-orbit split edges would be interesting. 
It is generally assumed that for isotropic light J° /^£ 2 = 
2:1. However, a strong polarization of the 5d electrons 
in the ground state would give deviations from the sta- 
tistical branching ratio. The 2p — > 5d circular dichro- 
ism certainly requires more experimental and theoreti- 
cal research. Explanation of circular dichoism by band 
structure calculations have so far only considered the df 
spin interaction. The difficulties in explaining the cir- 
cular dichroic branching ratios directly implies that the 
coupling of the 5d band to the local 4/ states is not well 



understood. Band structure calculations including or- 
bital polarization are necessary in the interpretation of 
the circular dichroism. 

The experimental developments in x-ray inelastic scat- 
tering are more recent. It is now rather well established 
that on-resonance the excitation and decay cannot be 
decoupled. For a delocalized system this implies that 
one has to take into account the crystal momentum of 
the core hole (although core-valence interactions in the 
intermediate state might change its value). For a lo- 
calized system one has to consider the angular momen- 
tum. Within an independent electron model the angu- 
lar momentum is conserved. However, interactions with 
the valence shell change the values of j and m and in 
general for an intermediate-state eigenstate one has to 
consider expectation values of, e.g., j z of the core hole. 
This quantity could be obtained by, e.g., measuring the 
difference between left and right circular polarization 
of the outgoing light. Unfortunately these experiments 
are rather complex. A determination of the quadrupo- 
lar moment of the core hole polarization seems more 
promising since it involves a measurement at two differ- 
ent detection angles. The presence of core-hole polariz^ 
tion has been demonstrated in resonant photoemissionL3 
(/" -> d n+1 -> c' 2 l n+1 E k , where E k denotes a photo- 
electron) which is formally very similar to the scattering 
process l n — > cl n+1 — ► c'l n+1 . Recent resonant Raman 
experiments on Co show the presence of higher moments 
in the core hole distribution. EZJ 

The x-ray scattering process l n — > cl n+1 — > l n en- 
ables one to study valence band excitations. Although 
for electric multipole transitions AS = 0, these excita- 
tions also include spin flips, since the spin is not a good 
quantum number in the intermediate state as a result of 
the large core-hole spin-orbit coupling. Recently, octet- 
sextet transitions have been observed in Gd 3+ (ground 

state s Sr)Si 

2 ' 

A relatively simple way to study resonant inelastic x- 
ray scattering is fluorescence yield, since it does not in- 
volve the detection of the energy of the outgoing photon. 
The spectra for fluorescence yield already show signifi- 
cant deviations from the XAS cross section. 



APPENDIX A: 

The tensor products in this paper are done with 3j- 
symbols which are up to a factor equivalent to those with 
Clebsch-Gordan coefficients 

[a',bn? = E^v(-D'- A - A '(_ Z A \ J) 

AX' V s / 
= {-l) l '[x]-^Y,Abx>C?t vx , ■ (Al) 

XX' 

When more than one superscript is present the last one 
gives the rank of the tensor; a tensor without superscript 
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has rank one. Furthermore, we often make use of the 
numerical factors 



ni x = 



(21) 



I I J y /(2l-x)\(2l + l + x)\ 



and 



n X yz — 



x y z 




(A2) 



(A3) 



The latter factor is zero for odd x + y+z. In certain cases 
it is convenient to have a "generalized" form for n xyz : 



= i 9 



{g-2x)\(g-2y)\(g- 2z)\ 
(<? + !)! 



I\ 2 



{g-2x)\\{g-2y)\\{g-2z)\\ 



(A4) 



with g = x + y + z. Note that n xyz is equal to n xyz 
for even g. The purpose of these factors is to remove the 
square roots from the expressions. These square roots are 
a result of the normalization of the 3j-symbols which is 
inconvenient when dealing with physical quantities, such 
as, operators and combinations of spectra. 

Some useful tensor products are: the inner product, 
[a 1 1 b l ]^ ) nJ Q 1 —a l -h l , the outer product, [a, b] 1 =--^axb, 
and the coupling of two spherical harmonics of the 
same vector [C'(k), C (k)] x =C x (k)(-) W* . In 
this paper we also apply the orthogonality relation 
(a ! -b')(c'' ■d l ')='£ x [x}[a l ,c l '] x -[b l ,d l '] x with [x]=2x+l. 
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